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I. INTRODUCTION 

The field of molecular electronics has greatly benefited 
from the possibility of formulating transport problems in 
the weak-coupling regime as rate equations for sequen- 
tial tunneling processes?^ Given the spectrum of the 
molecule and the tunneling-matrix elements, the differ- 
ent features of the steady-state current-voltage charac- 
teristics can immediately be mapped onto the energetic 
availability or non- availability of certain jump processes. 
Rate equations are the Markovian kinetic equations for 
the lowest-order expansion of the von Neumann equa- 
tion with respect to the tunneling Hamiltonian while ne- 
glecting any off-diagonal elements of the reduced density 
matrix. Such an approximation is well justified for non- 
degenerate systems, when the differences of the molec- 
ular eigenenergies are much larger than any tunneling- 
induced level shift or broadening. As molecules often 
feature geometric and thus orbital symmetries, the sys- 
tem Hamiltonian of a single- molecule junction shows de- 
generate levels. Due to the quantum-mechanical nature 
of the tunnel junction, the rate-equation description is in 
general inadequate for these systems. 

The problem of electronic transport through quantum 
nanostructures with degenerate levels has already been 
given some attention in different parts of the literature. 
The comprehensive review of Markovian master equa- 
tions by TimrrH^ shows the equivalence of different meth- 
ods and approaches used to derive master equations for 
weak-coupling problems. In Ref s . ITT1 and [T2| the coupling 
of a spin-degenerate quantum dot to ferromagnetic leads 
causes coherent dynamics described by the full master 
equation which significantly differ from the spectroscopic 
picture found in rate-equation treatments. The problem 
of using the rate-equation formalism for molecules with 
orbital symmetries has already been addressed in our pre- 
vious study on Jahn-Teller molecules.^ An implementa- 
tion of the full master-equation formalism for a genuine 
molecular-electronics problem is discussed in Refs.RHland 
1151 however, only for a very special and rather complex 
type of molecule. 

The dynamics of a master equation®! is, in contrast to 
a rate equation, much less intuitive as the tunneling elec- 



trons are allowed to jump into and out of linear superpo- 
sitions of the degenerate molecular states. The coupling 
to the continuum of states in the electronic reservoirs via 
virtual transitions generates an intrinsic dynamics on the 
molecule that is not related to real changes of the number 
of electrons in the system. The image of the dynamics 
as a succession of well-defined quantum jumps between 
the leads and the molecule, which renders the rate equa- 
tion so simple in its use and tempting for application, is 
declared void by the quantum mechanical nature of the 
degenerate system. 

The basic phenomenology of quantum transport 
through nanostructures with orbital degeneracies in the 
absence of vibrations and in particular the proper deriva- 
tion of a Markovian master equation for the treatment of 
near-degeneracies has already been investigated by us in 
Ref. [TSJ There we establish the "decoupling paradigm" , 
which states that for generic tunnel amplitudes, there 
is always a basis of the molecular Hilbert space where 
one of the (near-)degenerate levels is decoupled from the 
drain electrode. Below the double-charging threshold, 
this level is rendered a dark state in which charge is accu- 
mulated, and electronic transport across the nanostruc- 
ture is strongly inhibited. This coherent current-blockade 
is only partially lifted due to the tunneling-induced renor- 
malization of the isolated structure's levels, as by tran- 
sitions via virtual intermediate states in the source elec- 
trode, the electron is moved out of the dark state and 
allowed to tunnel to the drain electrode. A particularly 
appealing interpretation of this dynamics uses the picture 
of a pseudo-spin and the tunneling-induced renormaliza- 
tion as a pseudomagnetic field acting on that pseudo- 
spin P21 

The purpose of this article is the extension and dis- 
cussion of the model to single-molecule devices, mainly 
the relation of the physics caused by the orbital degen- 
eracies to the vibronic dynamics of the molecular cage. 
After having defined the most general model of a lin- 
ear electron-phonon coupling, we shall show that for the 
study of degenerate systems this can be reduced to two 
genuinely different types of molecular models, which sig- 
nificantly differ in their transport properties. One of 
them, which we shall term "Anderson-Holstein model" , 
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will simply show a superposition of a vibronic sideband 
structure and the already known coherent current block- 
ade generic for degenerate electronic systems. The ratio 
of charging energy to vibronic energy will be shown to 
characterize the steady-state current -voltage profile: for 
large charging energy, the vibronic sidebands will appear 
as peaks instead of steps and thus render the appearance 
of negative differential conductance a generic property of 
single-molecule junctions with degenerate orbitals. The 
other, a Jahn-Teller active model, will show a strong de- 
pendence of the electronic transport properties on the 
coupling to the leads. 

In the last section, we shall extend the basic model 
incorporate elements that arc related to possible experi- 
mental issues including the modifications of the transport 
properties due to slight breaking of the orbital degen- 
eracy, general linear electron-phonon coupling, and the 
presence of many modes in the electronic reservoirs. 

Not surprisingly, our results will reproduce certain ef- 
fects and feat ures, which have already been reported by 
several groups,^ 1 * 14 * 15 * 17 ! what we, however, do want to 
show is that by using a bottom-up approach and adding 
complexity to the models in several steps, we succeed in 
tracing the fundamental and generic physics back to the 
intrinsic properties of the master equation and are thus 
in a much better position to actually apply the theory to 
experiments. 



II. GENERAL PROPERTIES OF THE 
HAMILTONIAN 



The model of a single-molecule junction, as we con- 
sider it in this article, consists of three parts: the 
single molecule itself, H mo \, the source and the drain 
electrodes through which electrons are injected and ex- 
tracted, -f/ieads; and a tunnel-coupling between the two, 
Ht- The electrodes, a being an index for left and right, 
are modeled as spinless, non-interacting Fermi gases, 
ffioads = J2ka £ kci a c ka , in the wide-band limit, that 
is with constant density-of-states vq. The molecule, in 
order to distinguish it from quantum dot systems, is a 
discre te ele ctronic system coupled to a single vibrational 
mo de EELis] rj-^g e i ec ^ romc p ar t m our model is rather 

simple. We consider two degenerate levels that can be 
detuned from their energy £- by the gate-voltage of the 
molecular junction, a Coulomb interaction of strength 
U between the two, but no intra-molecular tunneling. 
A small energy difference between the two levels of the 
order of the tunneling-induced level shift can easily be 
included in the model and the thus resulting derivation 
of the master equation by application of the singular- 
coupling limitp^ modifies the equation only marginally; 
we shall return to this later in section HV Al The molec- 
ular vibrations are modeled as a harmonic oscillator of 
frequency to being coupled to the electronic degrees of 
freedom by -ffei-ph specified in the next section. The 



Hamiltonian of the molecule is thus 
-ffmoi = (£o + eV g )(n t +n;) + C/7i 1 .7i i +ftwfe t 6+i7 c i_ph. (1) 

The two degenerate electronic levels are, in analogy to the 
notation of spin, labeled |f) and In case we have to 
sum over the different levels, we switch to denoting the 
levels by \a) and its opposite \<j). The tunneling part, 
where from the beginning, we assume the amplitudes t aa 
to be independent of the electrons' wave vector, is 



(2) 



For notational convenience, we define a coupling tupel 
T := (r L t, Ti4, r Rt , r^) being derived from the respec- 
tive Golden-Rule expressions T aa := ^-vo\t aC r\ 2 - 



A. Electron— Phonon coupling and Degeneracies 

The general form of a linear coupling of the oscillator's 
coordinate to the charge number of a degenerate two- 
orbital molecule is 



H, 



cl— ph 



E 



(3) 



When the expression for i? e i_ p h is diagonalized in the 
space of degenerate electronic levels, which is necessary 
for the derivation of the master equation, the only non- 
trivial couplings are those to the excess charge and the 
charge difference, respectively, 

ffel-ph - XH(b^ + 6)(A (n t + nj + X z (n t - n^). (4) 

The canonical transformation generally used to elim- 
inate the electron-phonon coupling — the polaron 
transformation^ — induces a renormalization of the 
electronic eigenenergies — the polaron shift — which is 
proportional to the square of the electron-phonon 
coupling strength. The energy of the state |f) will thus 
be renormalized by (Ao + \ z ) 2 hui and the one of \l) 
by (Ao — \ z ) 2 huj, accordingly. The Hamiltonian in the 
polaron picture is, however, the starting point for the 
perturbative analysis of the von Neumann equation, 
and the question, whether we deal with a degenerate, 
a near-degenerate, or a non-degenerate system refers 
to this picture and the renormalized energies thereof. 
In case the electron-phonon coupling is assumed to be 
sufficiently strong, such that the polaron shift is much 
larger than the tunneling-induced broadening T, the 
system is effectively non-degenerate, and the canonical 
rate-equation formalism can be applied. If the polaron 
shift is of the order T, the electronic term has to be 
treated in the singular-coupling limit which will be 



touched briefly in section IV A For the setting with 
strictly degenerate orbitals that we wish to discuss 
here, the polaron shift of both levels has to be equal, 
which limits the possible choices of the electron-phonon 
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coupling in the above expression to either Ao or X z being 
zero. For strictly degenerate systems, we thus assume, 
without loss of generality, 

H e i- ph = \hu(tf +b){n^±ni). (5) 

Choosing the plus sign, we obtain a trivial generalization 
of the single-mode Anderson-Holstein Hamiltonianf^l 
This choice is therefore termed the "Anderson-Holstein 
molecule" . The minus sign, on the contrary, makes 
the electron -phon on coupling that of an E S3 b Jahn- 
Teller effect We call this model the "Jahn-Teller 
molecule". Certain transport properties of Jahn-Teller 
systems in the rate-equation regime are analyzed in 
Ref. H31 By absorbing the sign into orbital specific 
electron-phonon couplings X a , the Anderson-Holstein 
molecule is defined by — A4., whereas the Jahn-Teller 
case is given when A-j- = — Aj,. 

The polaron transformation renormalizes not only the 
single-particle energy £o ^ £o — \ 2 hu> but also the charg- 
ing energy U h-> U— 2 sgn(A-|-Aj,)?UJ. As in the Jahn-Teller 
case, both orbitals are shifted to opposite directions, oc- 
cupying such a molecule with two electrons will result in a 
zero net shift of the adiabatic potential, see Figure[TJ The 
renormalization of the charging energy U will therefore 
be positive; in contrast to the Anderson-Holstein model, 
the Jahn-Teller molecule induces a repulsive interaction 
between the two electrons in the polaron picture, not an 
attractive oneP^ 



B. Existence of Dark States 

As we already know from the treatment of purely elec- 
tronic structures™ the master equation can be under- 
stood best in a coordinate system of the degenerate lev- 
els, in which one of them is decoupled from the drain 
electrode. Then a dark state is formed and causes the sta- 
tionary current to be strongly suppressed. In the polaron 
picture, where the tunnel amplitudes are matrix valued, 
this property is modified. Consider first the case without 
electron-phonon coupling. We apply a unitary transfor- 
mation in the two-dimensional complex vector space "H c i 
spanned by the operators d^ and dj,, 

fdA _ ( cos9 e^sin6>\ / d T \ , . 

\d 2 ) '~ \-e~^sva6 cos(9 J [dj ' 

All parts of the Hamiltonian except the tunneling term 
are invariant under this transformation. The tunneling 
Hamiltonian becomes 

Ht=Y, c L ((*of cos 9 + sin 6)d 1 

ha 

+ (t a i cos 6 - t a1 -e~ lip sin 9)d 2 ) + h.c. (7) 

By choosing suitable angles 8 and (p, the second term 
of one of the above equation vanishes for at least one 
electrode, allowing the formation of a dark state: the 



basis of the coherent current blockade as we have ex- 
plained in Ref. I16I Under certain circumstances, that is 
for specific choices of the tunnel-couplings, the Hamilto- 
nian can be further simplified. If the condition t a ^cos9 — 
t a ^e~ lip sin 6 can be fulfilled for every a, the second term 
in Eq. ^ will vanish completely and the system will only 
couple the state di to the electrodes, thus reducing it to 
an effective single-level system with tunneling matrix el- 
ements 1 1 | 2 = \ta\ | 2 + |i Q ^| 2 . This condition reads 

ip = arg and 6 = arctan p?ii 5 (8) 

for all a. In the following, we shall assume real t aa for 
simplicity and accordingly set Lp = 0. 

1. Anderson-Holstein 

In the case of phonons, we have to distinguish the 
Anderson-Holstein and the Jahn-Teller case. In the po- 
laron picture where the tunnel amplitudes are matrix- 
valued t aa i — ^ t aa e~ x "( bl ~ b \ the second term of Eq. ^ 
reads for electrode a 

clj ai co S d(e-^-V -e-^~V)d 2 . (9) 

The electron-phonon coupling of the Anderson-Holstein 
molecule is characterized by the condition A^ = thus 
rendering the term ^ zero. In this model, the decou- 
pling paradigm and the emergence of a dark state works 
the same way as for purely electronic levels. In the de- 
coupling regime, such a molecule is equivalent to a single- 
level molecule. The Jahn-Teller molecule, however, has 
A-f- = — making a detailed discussion of the transition 
rates necessary. 

2. Jahn-Teller 

Transforming the molecular Hamiltonian into the po- 
laron picture essentially consists of shifting the adiabatic 
potential of the oscillator by \> r 2\ a £ osc . The Franck- 
Condon matrix element M$ q ' = \{q | e - A -( bt - b ) \ q ')\ 2 at 
row q' and column q itself is the overlap of the original vi- 
brational state \q) with the shifted oscillator's state \q'). 
The modulus of the matrix elements is independent of the 
direction of the shift, the sign, however, is not. A transi- 
tion between states that both have an even or odd num- 
ber of excited quanta still is independent of the direction. 
A transition from a state with an even number of quanta 
to one with an odd number of quanta, however, is sensi- 
tive to the direction. In the matrix e A(fct ~ b) - e~ x( - b ~ b \ 
only the matrix elements belonging to excitations of an 
even number of oscillator quanta cancel. In the pre-factor 
of di, the matrix t at cos 9e~ x ^- b) + t ai sin6>e A(i,t_f,) , 
however, there are no cancellations to be expected. Only 
in the case of equal coupling, t aa = 1 for all a and cr, 
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the matrix elements for excitations of an odd number 
of quanta cancel. In that special case, the Hamiltonian 
couples transitions with \q — q'\ even to <i 2 and transi- 
tions with \q — q'\ odd to d\. Since the both subsets 
of transistions change the oscillator's energy by either 
even or odd multiples of huj, only, the system decouples 
into two independent subsystems with disjoint spectra, 
which in the weak-coupling limit can be treated by rate 
equations. As simple and tempting a treatment by rate 
equations might seem, in this particular case, the rate 
equation would have two stationary solutions, such that 
the asymptotic dynamics would strongly depend on the 
initial state of the system at time t = 0. 

In the Jahn-Teller configuration, in contrast to 
the purely electronic model or the Anderson-Holstein 
molecule, no electronic level can be decoupled completely 
from a single electrode by a unitary transformation. As 
explicated above, only a certain subset of transitions 
can be decoupled, the others remain with finite transi- 
tion rates. Since this property renders the decoupling 
paradigm for the Jahn-Teller molecule void at first sight, 
it restricts the current-blockade to the voltage regime be- 
low the first vibronic sideband, but in return allows to 
evaluate the model even in the regime where Eq. [8] is 
satisfied and the Anderson-Holstein model splits up. 



III. TRANSPORT PROPERTIES 

A. General Phenomenology 

In the weak-coupling regime T <C k^T, Ku), most trans- 
port calculations using rate equations show results, where 
the Coulomb-blockade physics of the electronic levels is 
augmented by a vibronic structure and the appearance of 
vibronic sidebands in the current -volt age profile. 18 Sig- 
nificant diversions from this picture only occur in the 
case of very weak electron-phonon coupling,^ where 
the vibronic dynamics heat the molecule, or very strong 
coupling,^' where the huge displacement of the oscillators' 
adiabatic potential causes an exponential suppression of 



the low-bias current — the Franck-Condon blockade. 

In the presence of two degenerate electronic levels be- 
ing described by a master equation for the full reduced 
density matrix, we expect this paradigm to still hold true. 
The Coulomb-blockade physics of the electronic structure 
is enhanced by the coherent current blockade^ due to 
the presence of the dark state. Depending on the ratio 
U/huj, we shall observe vibronic sidebands in the form 
of peaks rather than steps on-top of the profile of the 
current blockade (U > hui) or the current blockade mod- 
ifying the first few vibronic sidebands (Huj > U). The 
numerical evaluation of the stationary current for the 
two generic molecular models shown in Figure [2] corrob- 
orates this reasoning. The Anderson-Holstein molecule 
complies well with this argument as the formation of the 
dark state is not influenced by the vibronic structure. 
By contrast, the inability of the Jahn-Teller molecule to 
completely decouple one electronic level shows up in the 
modification of the transport properties at the first vi- 
bronic sideband only. 

The role of the Lamb-shift contributions is in general 
the same as in the purely electronic case, yielding an ad- 
ditional intra-dot tunneling Hamiltonian of the form of 
a pseudo-magnetic field B, which connects the two de- 
generate states via virtual intermediate states in the elec- 
tronic reservoirs. The x-components of these fields in the 
pseudo-Bloch equation for S q are, due to the symmetry 
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(p - q)hu) + U + eV s - e 
1 - /(e-Mo) 



de 



(q - p)huj + eV g 



de 



(10) 



We assume the electronic bands of the leads to be wide 
enough to ensure Vj ^de ~ 0. All summands vanish 
simultaneously at eV g — fi a — hU, where, similarly to the 
purely electronic case, the current is completely blocked: 
1 = 0. 

Independently of U, there is a contribution of the vi- 
bronic excitations to the pseudo-magnetic fields, which 
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FIG. 2: (Color online): Stationary current as a function of the bias eV s d for various values of U /huj = 1, 2, 4, 8, 16, 32, 64, 128 — red 
to blue curves. The electron-phonon coupling strenght is A = 1.1, the tunnel coupling is T — (1, 1, 1, 2). (a) Anderson-Holstein 
molecule. For small U/fko, the coherent current-blockade is only seen before the first vibronic sideband. For higher values of 
U/fko, the sidebands appear as a peak-superstructure on the large-scale current-blockade. Since the current-blockade and the 
total inhibition of current flow at eV s d = U/2 is the dominating feature of this dynamics, the vibronic sidebands become peaks 
thus rendering the appearance of negative differential conductance a generic property of this model, (b) Jahn-Teller molecule. 
The coherent phenomenology in the stationary current is limited to the regime before the first vibronic sideband, as due to the 
a z type electron-phonon coupling, the decoupling cannot be performed for all vibronic excitations of a given electronic level. 



goes beyond the purely electronic model. A non- 
interacting molecule with U = 0, will still have finite 
pseudo-magnetic fields as the terms with p ^ q in Equa- 
tion ( fl0| ) do not cancel but yield a contribution that is 
sharply peaked at the vibronic resonances. This causes 
slight derivations of the line shape in the differential con- 
ductance from the rate-equation result -^hl oc -J^foc 
/(l — /), the derivative of the Fermi function. 



B. Strong Electron— Phonon Coupling 

Similarly to the rate-equation treatment, we expect the 
molecular models to show significant influence of the vi- 
bronic structure on the transport properties in the regime 
of strong electron-phonon coupling.^21 J n Figure |3j we 
show the stationary current of the Anderson-Holstein 
model at intermediate electron-phonon coupling A = 1.1 
in contrast to strong coupling A = 4, where the system 
is deep in the Franck-Condon blockade regime. 

In Figure [3] (a), we see the suppression of the low-bias 
current due to the Franck-Condon blockade. The addi- 
tional suppression of the stationary current at higher bias 
is not known from the rate-equation model and due to the 
intrinsic mechanism of the master equation with coher- 
ences. We have shown previously^ and also in this pa- 
per that by identifying an exact dark state in the model, 
we can find the stationary state of the system. If we 
can find a state that can be populated easily but leaks 
only very little probability, this state is still a candidate 
to acquire much population in steady-state and apt to 
dominate the transport physics. Applying this to the 
large-bias regime of the Anderson-Holstein molecule in 
the Franck-Condon blockade yields an explanation for 
the observed effect. Assume a large charging energy of 



several hto and bias and gate voltage chosen such that 
the singly occupied state is aligned slightly above the 
drain electrode and the doubly occupied state well below 
the Fermi level of the source electrode. Then the state 
| cr,0) is almost dark because of the exponentially sup- 
pressed Franck-Condon matrix elements for the transi- 
tion | <t,0) H> |0,0) at the drain. In the conventional 
rate-equation formalism this would not reduce the trans- 
parency of the device as one could occupy the second 
level and tunnel off to the drain by converting the large 
charging energy into vibronic excitations. In the coher- 
ent setting, however, we are allowed to decouple one elec- 
tronic state from one electrode, and although we have up 
to now only decoupled one state from the drain electrode, 
we could easily decouple one from the source electrode; 
let this state be \a). The Franck-Condon blockade thus 
obstructs tunneling to the drain electrode, and the decou- 
pling from the source electrode does not allow to doubly 
charge the device, hence reducing the conductance of the 
device drastically. As before, the Lamb-shift contribu- 
tions help to regularize this picture, by transferring the 
charge from |er) to \<r), allowing a second electron to oc- 
cupy the device and converting the charging energy to 
vibronic energy in the tunneling to the drain electrode. 
In this particular setting, the roles of source and drain 
electrode have been reversed in the game of decoupling 
and transferring charge via virtual intermediate state of 
the reservoirs compared to the low-bias regime. 



C. Particular Phenomenology of Jahn— Teller 
molecules 

The electron-phonon coupling of a Jahn-Teller 
molecule is proportional to cr z in the electronic Hilbert 
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FIG. 3: (Color online): Stationary current of the Anderson-Holstein molecule in the Franck-Condon blockade regime for A = 4 
(a) and for intermediate coupling A = 1.1 (b). Although the obvious effect of strong electron-phonon coupling is a suppression 
of the stationary current at the Coulomb-blockade line and the double-charging threshold, the current profile develops sharp 
spikes at larger bias. The parameters are U = Whoo and T = (1, 1, 1.5 2 , 0). 



space and thus breaks the rotational invariance of the 
molecular Hamiltonian, which we have exploited to un- 
derstand the steady-state dynamics. We have already 
shown that this results in the inability to decouple an 
electronic level completely from one of the electrodes, 
thus obstructing the use of generic tunnel couplings. 
Looking closer at the quantitative properties of the mas- 
ter equation, the consequences can be seen in a number 
of effects. 



The different directions of the shift of the adiabatic 
potentials implied by the polaron transformation intro- 
duce a sign dependence on certain rates in the master 
equation. In particular the rates "fP q , will have dif- 
ferent signs for \p — q\ being even or odd. When we 
encounter a situation, where the s-component of the 
pseudo-magnetic fields plays an important role in the dy- 
namics, for example by transferring charge from a dark 
state into a conducting one, the structure of the Jahn- 
Teller molecule's electron-phonon coupling will cause a 
significant reduction of B x and thus of the current restor- 
ing force. The consequence is, in comparison with the 
Anderson-Holstein molecule, a profile of the coherent 
current blockade, which hardly shows an influence of the 
Lamb shift at all, having a much higher contrast as is 
shown in Figure [4j 



A second effect of the broken symmetry in the Jahn- 
Teller molecule is the ability to push the system into 
the parameter regime, where the Anderson-Holstein 
molecule would naturally split into two independent sys- 
tems, that is, when the ratio t a ^/t a ^ is independent of the 
lead index a as given by Equation [8j The main difference 
between the two models is that Jahn-Teller molecules the 
pseudo-magnetic fields are non-zero. In this regime, the 
pseudo-Bloch equations for the model in a voltage regime 
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FIG. 4: Effects of small pseudo- magnetic field in the x- 
direction. Tunnel coupling T = (1, 1,0, 1), zero gate voltage, 
Anderson-Holstein: red, Jahn-Teller: blue, U = 2.35/iw 
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Because B x =^ 0, this equation is uniquely solved with 
5 = 0. A completely vanishing pseudo-spin is not only 
the statement that the coherences of the stationary den- 
sity matrix are evaluated zero but that also the steady- 
state populations of both electronic levels are equal. Such 
a result is unexpected insofar as it is independent of the 
ratio L^/r^, which in a rate-equation treatment would 
determine the pseudo-spin's z-component. Also, a van- 
ishing pseudo-spin means that the reduced density ma- 
trix in the singly-charged sector is proportional to the 
unit matrix and invariant under rotations in the elec- 
tronic Hilbert space. 
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IV. MODIFICATIONS OF THE BASIC MODEL 

In this section, we shall discuss several extensions of 
our basic model of a single-molecule junction with or- 
bital degeneracies, which are related to experimental is- 
sues such as not having exact degeneracy or single-mode 
reservoirs. 



A. Slight Breaking of the Degeneracy 

The requirement of strict degeneracy of the two or- 
bitals is quite unphysical, when one takes into account 
all the physics of a single- molecule junction that we have 
neglected in the abstract model. However, the deriva- 
tion of a master equation for the reduced density ma- 
trix requires strict degeneracy in order not to sacrifice 
the density matrix's coherences to the rotating-wave ap- 
proximation. Although there have been some attempts 
in the literature to derive a proper master equation for 
this case via th e, not necessarily positive, Bloch-Redfield 
equation) 11 * 15 ^ we have shown in a previous workP^ how 
such an equation can be derived rigorously by treating 
the breaking of degeneracy as a perturbation. 

Assume the degeneracy of the electronic system Hamil- 
tonian Hg — (eo+eV e ) d\d ai or using the pseudo-spin 
notation H — (e g + eV^Id, is broken by Hq :— ^ilcr z , 

that is Hq = jQ(dtd<f — rftdj.) in second quantization. 
Assuming this perturbation to be a slight breaking of 
the degeneracy that will not destroy the coherences in 
the reduced density matrix implies the relation $7 ~ T <C 
k^T, hu. Following the derivation given in Ref. 1161 we 
only have to add the Liouvillian Cn = [Hq, •] to the mas- 
ter equation found for the degenerate, f2 = 0, model. 
This term induces a second contribution to the Hamilto- 
nian part of the master equation besides the Lamb shift; 
the splitting essentially adds to the z-component of the 
pseudo-magnetic field, but in contrast to the Lamb shift, 
fHs a free parameter of the theory. Its influence on the 
steady-state transport properties of our single-molecule 
devices will be investigated in the next few sections. 



1. Anderson-Holstein Molecules 

As we have shown before, for intermediate electron- 
phonon coupling, the Anderson-Holstein molecule's be- 
havior is a superposition of the phenomenology of the 
electronic levels and the vibronic sideband structure. 

In the numerical evaluation of I(eV s< i) as a function of 
$7 in Figure [5] (a), the current suppression is lifted for 
larger splitting f7, and the small vibronic peaks develop 
into well-defined steps of the current profile, which are 
typical for the rate-equation treatment. Due to different 
tunncling-induced pseudo-magnetic fields B q , for each vi- 
bronic excitation but constant splitting $7, the lifting of 
the coherent structure and the recovery of the flat profile 



is different for each sideband. In Figure[5](a), the various 
vibronic sidebands emerge at different values of $7. 



2. Jahn-Teller Molecules 

Due to the inability to decouple all vibronic excita- 
tions of a single electronic level from the drain electrode 
in the Jahn-Teller model, the coherent current blockade 
being the generic phenomenology of Anderson-Holstein 
molecules is in general only visible if the zero of the 
source's pseudo-magnetic field is at a voltage below the 
first vibronic sideband. In Figure [5] (b), we show the 
fi-dependent stationary current at zero gate voltage for 
intermediate charging energy. As we have claimed, the 
current blockade is localized close to zero bias and van- 
ishes quickly as the splitting $7 is increased. Since the 
zero of Bg, where the current blockade is fully developed, 
is shifted towards higher voltages for large charging en- 
ergy and therefore far beyond the first vibronic sideband, 
systems with larger Coulomb repulsion will not show any 
suppression at all. 

For finite f7, however, we can shift the zero of i?s so 
far that it moves below the first vibronic sideband al- 
lowing the vibronic ground state to become dark. We 
show our findings in Figure [6j where we plot the cur- 
rent suppression due to finite $7. Figure [6] (a) shows the 
current profile for eV s = —hiu and varying splitting Q. 
The dip in the stationary current for 17 rs 0.6 is due to 
the mentioned effect. Diagrams (b) and (c) show a com- 
plete current-voltage profile in the (eV g , eV r s d)-plane. For 
$1 = (b), the current suppression is absent for all volt- 
ages. For 17 = 0.6 (c), on the contrary, there is a deep 
trough of suppressed current, slightly reminiscent of the 
phenomenology of the Anderson-Holstein molecule, but 
only for the first vibronic sideband and negative gate 
voltage. 



3. General linear electron-phonon coupling 

In the introduction, we have considered the general 
situation of linear electron-phonon couplings. For de- 
generate electronic systems, these can be always diago- 
nalized leaving only the identity and the a z component 
of the electron-phonon coupling matrix in the electronic 
Hilbert space 'H e \ A = Aold + X z a z . Due to the polaron 
transformation, the system's electronic levels will acquire 
an energy difference proportional to 4AoA 2 . The strictly 
degenerate theory, $7 = 0, can only be applied to either 
the Anderson-Holstein model, X z = 0, or the Jahn-Teller 
model, A = 0. Using the concept of near-degeneracies, 
however, we can go beyond this sharp distinction and 
consider molecules with both Ao and A z being non-zero 
and include the splitting due to the polaron transforma- 
tion as a near-degeneracy $7. Such models are interest- 
ing, because the different electron-phonon coupling of the 



FIG. 5: Scan of Q at eV g = 0, eV^d £ [—5fkjj,5huj], U = 2.35 huj, (a) Anderson-Holstein molecule with A = 1.1, T = 
(1, 1.5 2 , 1.5 2 , 1). The splitting is not measured in units of hu but, due to the singular-coupling limit Q, ~ T in units of ^f^o- 
(b) Jahn-Teller molecule for the same parameters. 
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FIG. 6: (a) Scan of fi for the Jahn-Teller molecule at eV e = -huj, F = (1.7 2 , 1, 1.4 2 , 1), U = lOfioj, A = 1.1, and knT = 0.02 fiw. 
(b) stationary current for f2 = 0, eT4d € [0, 5foj], eV^ G [— 2.5 hui, 2.5 hu)]. (c) same plot, however, with Q — 0.6, where similarly 
to the Anderson-Holstein molecule, a trough in the stationary current due to the coherent blockade effect is visible. 



system's levels amount to different Franck-Condon ma- 
trices, for example one level being weakly coupled, the 
other being in the Franck-Condon blockade regime. 



B. Multi-Mode Reservoirs 

In a number of recent publications, electronic trans- 
port through single molecules is realiz ed by using sus- 
pended carbon nanotube quantum dots P^l Following the 
experimental set-up used in Ref. |H1 a top-gate attached 
to the suspended nanotube defines the central piece of 
it as a quantum dot, the outer pieces as electron sup- 
plies directly connected to the metallic source and drain 
electrodes, and all three regions being separated through 
mechanical deformations of the nanotube. The princi- 
pal effect of the thus created tunnel barriers to the left 
and the right of the nanotube quantum dot is that we 
cannot, in general, assume the valley quantum number 



of graphene to be conserved in tunneling anymore. If 
we denote the valley quantum number of the quantum 
dot region in accordance with the notation for degener- 
ate orbitals used throughout this article by a and the 
valley quantum number of the supplying nanotube elec- 
trodes by r, the tunnel matrix elements acquire an addi- 
tional index Y a ar . Instead of formerly four, we now have 
eight tunnel matrix elements, and it is immediately clear 
from this number that the paradigm of electronic decou- 
pling at the drain electrode, which we have established in 
the discussion of the Anderson-Holstein molecule, is no 
longer generically true for such systems. The reduction 
of the master equation to a rate equation in a suitably 
chosen basis of the electronic Hilbert space can, how- 
ever, nonetheless be achieved. Let us consider the non- 
interacting retarded self-energy of the tunneling problem 
with respect to a single electrode, where we assume that 
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FIG. 7: (a) Two-mode reservoir interacting with a two- 
level quantum dot with U = 2.35 hui, knT — 0.05 hu). (b) 
Two-mode reservoir interacting with an Anderson-Holstein 
molecule. U — 10 Hcj, knT — 0.05 htu, and Af = Aj, = 1.1. 
The coupling to the left electrodes are equal and set to 
unity, and the ones to the right electrode are £r-|-,i- = 1-4, 
tat, J. — £iu,T = 0- The coupling tru^ is scanned from to 1.4 
(red to blue curves) in order to increase the eigenvalues of £r 
from (£~,E + ) = (0,1.4), where the effect of the two-mode 
reservoir is absent, linearly to (E~,E + ) = (1.4,1.4), where 
the coherent current blockade is gone. 



the electrodes' distribution function is diagonal in r, 



oc 



X/r X/t 7a, t 

X^t 7a, t X/r ^ar 



(11) 



This matrix is hermitian and can thus be diagonalized 
with its off-diagonals being zero. If we choose a to be 
the drain electrode, the master equation without the 
pseudo-magnetic fields is rendered a rate equation. But 
now there is no generic dark state and transport does 
not proceed via just one electronic levels but including 
both; however, in a manner known from rate-equation 
theory. Assume the eigenvalues of Ed are E + ^> E~. 
Then the electronic system will accumulate more popu- 
lation in the eigenstate of E _ and its exit rate towards the 
drain electrode will determine the steady-state current. 
The pseudo-magnetic fields induce an intra-dot tunneling 
term, which pushes population from ||) to |f) and thus 
increases the steady-state current when the bias voltage 
is tuned away from the point \B\ = 0. In Fig.[7j we show 
numerical results for this model, where we increase the 
smaller eigenvalue of Ed from zero, which corresponds 
to the situation found in single- mode reservoirs, to E + , 
where the negative differential conductance structure has 
vanished, for in that case, both reservoir modes are cou- 
pled to the molecule equally well, and the drain's self- 
energy is proportional to the unit matrix. 

The physical significance of this result is the following. 
The graphs in Fig. [7] (b) do resemble the current -voltage 
diagrams computed in Ref. 1111 And indeed, the problem 
of electronic transport with ferromagnetic leads formally 
resembles our model as there it is the real electronic spin, 
which is present both in the electrodes and the quantum 
dot, but due to the ferromagnetic ordering not conserved 
in the tunneling. The primary influence of the presence 
of more than one electronic mode in the leads is the lift- 
ing of the strict current-blockade and the accompanying 
softening of the negative differential conductance at the 
vibronic sidebands. 



V. CONCLUSION 

The transition from the rate-equation physics of a sin- 
gle spin-degenerate quantum dot to a single-molecule 
junction is essentially the addition of inelastic transi- 
tions, which first of all cause vibronic sidebands. Only 
in certain regimes of the model's parameter space, ad- 
ditional vibronic physics like for instance the Franck- 
Condon blockade are visible. In the present article, 
we have proceeded in a similar manner to extend the 
master-equation theory for systems with orbital degen- 
eracies from electronic levels to molecular models, where 
the main persistent phenomenology is the addition of vi- 
bronic sidebands to the generic current suppression due 
to the coherent interaction of degenerate or near degener- 
ate electronic levels. With two instead of only one level, 
we have shown that already in the Hamiltonian, the de- 
coupling paradigm, being the essential tool to understand 
the transport dynamics, is modified by the matrix struc- 
ture of the tunnel amplitudes in the polaron picture. We 
have shown that in the strictly degenerate case, there are 
two generic models, which we have termend Anderson- 
Holstein and Jahn-Teller, with opposite phenomenology. 
While the Anderson-Holstein model complies well with 
the idea that an electron-phonon coupling simply adds 
the well-known vibronic physics to the steady-state dy- 
namics of the electronic model, we have shown that the 
phenomenology of the Jahn-Teller model shows many 
different transport regimes. The main reason being the 
inability to completely decouple an electronic level from 
one electrode by unitary transformations of the degener- 
ate orbitals. The decoupling can only be achieved for a 
certain subset of vibronic excitations, which significantly 
influences the steady-state properties of the transport 
model. 

We have thoroughly discussed the phenomenology of 
the generic models in various parameter regimes both for 
the tunneling amplitudes and the electron-phonon cou- 
pling. In the last part of the paper, we have relaxed the 
requirement of strict degeneracy of the electronic orbitals 
and allowed for so-called near-degeneracies, whose influ- 
ence on the transport properties of purely electronic sys- 
tems we have already investigated in a previous publica- 
tion. We have applied our results to the vibronic models 
and shown how electronic transport of both models was 
modified in these regimes. As a last step, we have turned 
to the electrodes and discussed the effect of multiple elec- 
tronic modes in the leads. Our main result is a modi- 
fication of the decoupling paradigm used in our theory. 
Although the coherent current blockade being generic for 
our models is lifted and, depending on the parameters, 
only slightly reduced by the presence of more than one 
electronic mode in the electrodes, applying unitary trans- 
formations in the space of the degenerate orbitals of the 
molecule again proves a suitable tool to reduce the mas- 
ter equation to a rate equation for specific voltages and 
thus understand its stationary solution intuitively 

We are well aware that none of the models, we have 
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discussed so far, would be a complete and quantitative 
description for any experiment as the models are too ab- 
stract. However, we have shown that already our models 
show the generic phenomenology found for more com- 
plex systems, for example in Ref. 1151 and that we can 
understand its physics on a very fundamental basis. We 
have shown how to bring the master equation into a form 
where its dynamics can be understood intuitively using 
basic principles already known from rate-equation theory. 
With our approach to the theory of electronic transport 
through single-molecule devices with orbital degeneracies 
we therefore provide a basis for future research in the 
physics of coherent interactions in sequential tunneling 
beyond the rate-equation approach. 
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Appendix A: Derivation of the Master Equation 

Our aim is the study of the steady-state transport 
properties of the molecular models in the weak-coupling 
limit, which amounts to only consider sequential tunnel- 



ing processes. To achieve this, we will derive a marko- 
vian master equation closely following RefsP^S This 
derivation, whose principal ideas have already been used 
by us in Ref. 1161 is not contained in the extensive review 
by Timmp^ but in our opinion it is the most intuitive 
and rigorous derivation, because it reduces the necessary 
assumptions by a large degree when compared to other 
approaches. The Hamiltonian of our problem is that of 
a system-bath interaction and has the general form 

Ht = H s +H E + £ff S -E. (Al) 

There is a finite-dimensional discrete system H$, the sin- 
gle molecule, an infinite particle reservoir, the environ- 
ment -ff E , which in our model are the electronic leads, 
and a system-bath coupling -ffs-E; namely the tunneling 
Hamiltonian. The parameter £ in front of i?s-E symbol- 
izes the weak-coupling assumption. Later £ — > will be 
implied and, in order to obtain a non-trivial result, an 
asymptotic time scale will be chosen. The quantity of 
interest is the reduced density matrix of the system S, 
ps := Tr E (p). obtained by means of a projection operator 
Vp := Tr E (p)®PE = PS®Pv- Also Q := Id-V. The en- 
vironment's density matrix will be the equilibrium distri- 
bution, which for fermions is just the Fermi-distribution. 
The von Neumann equation p = —t\H^,p\ =: —iC^p is 
split into an equation for Vp and one for Qp and in- 
tegrated formally by using the variation-of-constants for- 
mula known from the theory of ordinary differential equa- 
tions, where Ci refers to the Liouvillian with respect to 
Hi, 



ps(t) = e~ l£s Vs(0) - e J e~ lCs ^ lj' Tr E (Cs^e^^ £ S - E p E ) Ps(u)duj ds. (A2) 

By choosing a factorized density matrix as the initial condition p(0) = Vp(0), we implement the Born approximation 
and obtain the integral equation for the reduced density matrix alone. The substitution s — u + v allows to pull 
the reduced density matrix out of the inner integral, change the bounds of integration, and move to the interaction 
picture with respect to the system Hamiltonian 

Ps (*) = Ps(0) ~ ( 2 J q e~ lCss ' e* Csv Tr E (z^'^/We) dv\ e^p^s) ds. (A3) 

We aim at deriving a Markovian master equation for which we have to consider a time scale on which all memory 
effects are absent. This time scale is defined by the asymptotic time r = £ 2 i, which is being held constant in the 
limiting process £ — > 0, such that 

Ps(r) = Pl(0) - i 2 £ e Csfl e* Csv Tr E (£ S - E e-^£ s _ E p B ) dv J e"^ pl(s) da. (A4) 



Proceeding with the limit £ — > 0, we find that the oper- An expansion of the Liouvillians in terms of commutators 
ator in curly brackets converges^ to yields the nested commutator structure well known from 

standard perturbative treatments of the problem. 10 The 



V:= y°°e i£sI 'Tr E (c^e-< Cs+c ^ v £ s _ E p E ) dv, 

(A5) 
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action of the two exponential factors exp(±z£sc/£ 2 ) to 
the left and right of T> is to produce the time average 
T> of the operator.^ The Markovian master equation is 
thus 



Ps(t) = -T>p s (t), 



(A6) 



whose explicit representation for the models discussed in 
this paper is given in Appendix |B 1| In Refs. [55] and 
1271 one finds a very elegant argument, why the time av- 
erage of T>, being an exact result in the limit £ — > 0, 
yields the secular or rotating-wave approximation, which 
is usually applied to decouple the coherences between 
non-degenerate states from the respective populations in 
the reduced density matrix. 

We compute the stationary current within the same 
framework and start from the general definition of the 
observable. The stationary current through lead a is 



that cancels the first term, and changing to the asymp- 
totic Markovian time scale, which shifts the integration 
boundaries to infinity and takes the density matrix out 
of the integral. The Liouvillian VCs+eQ acting on a fac- 
torized density matrix is ££s-e- The emerging formula is 
identical to the standard, perturbative derivation, where 
the factorization of the density matrix has to be put in by 
hand but only after certain manipulations have been ap- 
plied to the formula. However, we now know that we can 
use the already obtained density matrix from Eq. (A6) 



and put it into the formula for the stationary current. An 
explicit formula in the basis of the single-molecule device 
is given in Appendix |B 1| 



Appendix B: Explicit Form of the Master-equation 



\Ia) — Trs +E (/ a ps+E)- 



(A7) 



The quantum mechanical current operator for lead a 
is the time-derivative of the number operator of this 
electrodeP^ 



I a = 



ha 



(A8) 



which is the 



the real part of this operator. We split the expectation 
value (I a ) into two parts by inserting Id = V + Q, 

Tr s+E (/ Q Pp s+E ) + Tr s+E (/ Q Qp s+E ). (A9) 



imaginary part of the operator 
The tunneling Hamiltonian Ht is 



a, 



The first term naturally vanishes, because V p = (Tr E p)® 
p E , and the system-bath coupling has zero expectation 
value with respect to the equilibrium distribution of the 
electrodes. It is the projection onto the complement 
which is the interesting term. The time evolution of 
Qps+e is given by 23 



^-iQC s+E Qt 



Qps+e(0) 



> f d S e- lQCs + EQ ^QC s+E Vp s+E ( S )- 
Jo 



After integrating the defining equation for the current 
I = —N being the time derivative of the number operator 
of the respective electrode, we obtain an integral equa- 
tion for the latter. We apply the same manipulations 
as before, using a factorized initial condition for ps+E 



1. Master equation 



In this section, we give the explicit representation of 



the master equation ( A6 ) for a two-level molecule with 
arbitrary Coulomb interaction U. The symbols used 
are defined by p q := |0, g) {Q,q\, p% := \a,q) (a,q\, and 
v\ '■= It-Jo l) CN-> <l\ > where the first entry specifies the 
charge state of the molecule and the second the number 
of excited quanta of the harmonic oscillator in the po- 
laron picture. The electrons are fermions, which requires 
to choose a definition of the wavefunction of the dou- 
bly occupied state, |2) := = d|d||0). The matrix 
elements of the molecular terms in the tunneling Hamil- 
tonian therefore differ for having a neutral or a doubly 
occupied state. 



i o) 



7t 



I) — tal — (t I Ct"t 



We also assume only real tunnel matrix elements and 
a wide-band limit of the electronic reservoirs, that is 
the density of states is energy independent and then 
set to unity for convenience. We define several short- 
hands: t^ a is the tunnel amplitude in the polaron picture 

t a a (q | e~ x °( h | p), with A<j incorporating the direction 
of the adiabatic potential's shift of the specific model. 
1^ := lei 2 , 1™ := t%t% = t%t%. A missing index 
on r or 7 indicates that the sum over the respective index 
is implied. The symbol Vj indicates the principal value 
integral. 



pg = - f^ p i 4)Wo + 2?r - Uz p i - 4)) ( r K 



a, p 



7S 9 2RepS 



(Bl) 



12 



P% 



2^(1 - f a {el - e\))Y qp v \ + 2n ]T f a {e\ - e\){YZpl ~ if 2Re^ ff 



(B2) 



2^£(iU £ 2 " e?)rSI + (1 - /«(ef - el))TZ)pl - 2 7 r£(-/ Q ( £ f - e\)Z q + (1 - / a (e? - eg))^) Bep* s 



2^E(/«( £ ? ~ £ o)« + (i - /a(4 - ^))rM - 2P /efe 



/a(E) 



1 ,,,, I ■ 



S-t En £ 



(B3) 



pl s =2nJ2(Usl - eg) 7 « - (1 - f a (4 ~ 

- * E [(-/«( £ 2 - e?)7S 9 + (1 - U{e\ - el))^){p% + p%) + (f a (4 - e\)T^ + (1 - f a {e\ - 4WS)pU 

- <p fa{ J a{ $ A (Kl - t%)pU + i p a q (pl - pD) + J"¥ e) (TO - r^) P L + t?(pI - 

(B4) 



Some cosmetic simplifications can be applied by noting 
that T pq = T qp and abbreviating the Fermi factors f pq := 
f( £ i - 4 - Pa) and fH := f(e p 2 -e\- fi a ) = f{e\ + U- 
Eq — fj, a ). In the Anderson-Holstein case, \ a = X s , also 
ry-pq = ryip^ The stationary current through lead a can 
either be computed using the formulae given in section [A] 
or derived from the above equations of motion by noting 
that this current is given by the a-contribution of the 
expression J2 P (Po ~ Pf)- Either way we find 



27T 



^(ia) = E f'lWo -EC 1 " ^)rK 

-E[( 1 -/S) r --/Mp? 

-E^-ZM+ZMSRe^. (B5) 
2. Pseudo-Bloch representation 



The master equation ( A6 ) for degenerate two-level sys- 
tems can be cast into a more intuitive form, since the 
electronic system admits the description by a pseudo- 
spin 5 for the singly charged part of the reduced density 
matrix. Due to the vibronic structure there is a pseudo- 
spin for every excited vibrational state |<j; q) there is an 
additional index S q . The components of the pseudo-spin 
are naturally defined by 



Due to the possibility of changing the charge state, the 
modulus of this spin is not constant, and the trace of 
Pi Pi '■= Pcr<? as t ne zeroth component of the Pauli- 
matrix representation has to also to be included. Here, 
we give an explicit representation of the master equation 
in terms of the pseudo-spin S q and the populations of the 
three charging states. 

& =tt £ + (/» - (1 - fZ)n q M 

a.p 

2(1 - Or«] - £ tt XI + (1 - /£)] T q J>S q 

a 

- (B pq + Bf) x S q . (B7) 

The pseudo-magnetization n vq is defined by the tunnel 
couplings, 



ffPQ 



(2 7 w o,ir-r§« 



(B8) 



The pseudo magnetic fields are defined by the principal 
value terms 

Sf:=E^/ - «tefg, (B9) 

a J b i fc i fc 

B?:=Y. V i y f / £) den q J>. (BIO) 



S q := (Rep^,Imp q n ,p q n - p q u ). (B6) 



The equations for the populations are in matrix form 
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By considering the dynamics of the electronic levels only 
and neglecting any oscillator excitations, the summations 



over p cancels due to 5 Pyq and we obtain equations similar 
in form to Ref. ITT1 
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